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ABSTRACT 


A  material  model  for  soft,  wet  soil  was  generated  to  simulate  the  deformation  behavior  of  thawing  soil 
under  vehicle  loading  on  paved  and  unpaved  roads.  Freeze-thaw  action  produces  a  loose,  wet  soil  that 
undergoes  large  deformation  when  subjected  to  vehicle  loads.  The  soil  modeled  is  a  frost-susceptible  fine 
sand,  which  was  used  in  full-scale  tests  of  paved  and  unpaved  road  sections  in  CRREL’s  Frost  Effects 
Research  Facility  (FERF).  The  soil  was  subjected  to  a  full  suite  of  saturated  and  unsaturated  triaxial  test¬ 
ing,  using  density,  moisture,  and  loading  conditions  duplicating  those  experienced  during  the  freeze-thaw 
testing  in  the  FERF.  Material  parameters  were  generated  for  a  capped  Drucker-Prager  plasticity  model. 
These  were  calibrated  in  triaxial  test  simulations  using  the  commercial  finite  element  code  ABAQUS. 
The  material  model  was  then  implemented  in  several  three-dimensional  finite  element  simulations  for 
validation  and  robustness.  The  model  for  Lebanon  Sand  was  compared  to  the  same  model  for  other  granu¬ 
lar  materials. 
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1  INTRODUCTION 

The  Problem  of  Spring  Thaw 

Spring  thaw  is  a  critical  time  of  year  for  the  deterioration  of  roads  and  air¬ 
fields.  The  prediction  of  performance,  or  deterioration,  of  pavement  structures  is 
important  not  only  for  maintenance  scheduling,  but  also,  for  the  military,  for 
deployment  and  operations  at  all  times  of  year.  In  developing  nations  the  use  of 
secondary  roads  and  airfields  for  rapid  deployment  will  be  critical,  and  the  per¬ 
formance  of  these  structures  must  be  optimized. 

The  types  of  deterioration  or  distresses  most  common  during  spring  on  sec¬ 
ondary,  unpaved  roads  are  rutting  and  roughness  (Shoop  et  al.  2002),  as  illus¬ 
trated  in  Figure  1.  Deformation  of  the  road  surface  (either  paved  or  unpaved) 
during  spring  thaw  is  nearly  always  the  result  of  the  deformation  of  the  thawing 
layer,  which  has  undergone  a  reduction  in  density  and  is  often  very  wet  or  satu¬ 
rated  as  a  result  of  the  freezing  process.  The  deformation  of  the  loose,  weak, 
thawing  layer  is  largely  plastic,  consisting  of  both  compaction  and  shear. 

Because  the  thickness  of  the  thawing  layer  varies  throughout  the  thawing  season 
(  during  spring  or  intermittent  thaws),  a  method  to  simulate  the  impact  of  both  the 
plastic  deformation  and  the  varying  layer  geometry  was  desired.  Therefore,  we 
developed  a  finite  element  model  of  the  layered  thawing  road  system.  The  pri¬ 
mary  challenge  was  to  capture  the  deformation  behavior  of  the  thawing  layer, 
which  had  not  been  previously  documented. 
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Figure  1.  Rutting  on  unpaved  secondary  roads  during  spring  thaw.  (Photo 
by  George  Blaisdeii,  Hanover,  NH.) 


Approach  to  Predicting  the  Effects  of  Freeze-Thaw  on  Vehicle  and 
Pavement  Performance 

Because  of  the  difficult  and  time-consuming  nature  of  constructing  full-scale 
test  sections  and  field  experiments  on  thawing  ground,  we  aimed  to  create  a  finite 
element  modeling  capability,  validated  with  experimental  data,  that  would  then 
be  useful  for  performing  computer  experiments  on  a  wide  range  of  off-road, 
pavement,  and  vehicle  or  aircraft  loading  situations.  To  do  this,  we  first  had  to 
develop  a  material  model  to  represent  thawing  soil,  as  none  was  available. 

This  paper  documents  the  experimental  development  of  a  constitutive  model 
representative  of  a  typical  frost-susceptible,  thaw-weakening  material.  The  soil  is 
well  characterized,  with  additional  data  provided  in  the  appendices.  Although 
other  material  models  were  tried,  the  material  model  chosen  was  a  modified 
Capped  Drucker-Prager  constitutive  model  as  implemented  in  the  commercial 
finite  element  code  ABAQUS.  An  extensive  laboratory  program  to  fully  charac¬ 
terize  the  mechanical  behavior  of  the  soil  was  completed.  The  test  methodology 
and  results  are  presented,  along  with  the  material  model  and  the  calibration  of  the 
model  parameters.  Finally,  the  thawing  soil  material  model  is  successfully 
implemented  in  a  3D  finite  element  simulation  of  a  wheel  rolling  on  an  off-road 
soil  and  a  layered  pavement  using  a  model  incoiporating  freeze-thaw  layering. 
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2  METHODS 


Material  Model  for  Thawing  Soil 

The  deformation  of  the  layered  pavement  system  during  spring  thaw  is 
largely  (perhaps  entirely)  the  result  of  deformation  of  the  weak,  thawing  layer. 
The  thawing  soil  may  lose  compaction  as  a  result  of  excess  water  being  drawn  to 
the  freezing  front,  expansion  during  freezing,  and  ice  lens  formation.  Also,  the 
loose  material  is  often  wet  or  saturated  during  thaw  when  water  is  trapped  in  the 
thawing  layer  by  the  impermeable  frozen  layer  below.  This  material  deforms 
plastically  in  both  compaction  and  shear  (Shoop  et  al.  2003). 

The  major  deformation  modes  of  the  thawing  layer  must  be  adequately  repre¬ 
sented  by  the  material  constitutive  model.  A  Capped  Drucker-Prager  (CDP) 
plasticity  model  was  chosen  to  capture  the  elastic  and  plastic  behavior  of  the 
material  under  vehicle  loading  during  spring  thawing  conditions  (i.e.  the  soil  may 
be  loose,  the  water  content  is  wet  of  optimum  to  fully  saturated,  and  the  loading 
is  short  term).  The  features  of  the  CDP  model  include  elastic  response, 
hardening/softening  as  a  function  of  volumetric  plastic  strain,  inelastic  behavior, 
and  yield  behavior.  In  the  CDP  model  as  implemented  in  the  commercial  finite 
element  software  ABAQUS  (HKS  2002),  the  shear  line  is  fixed  ( d  and  (i  are  con¬ 
stant)  and  hardening  (  compression  of  the  material)  occurs  on  the  cap.  The 
parameters  include 


d 

P 

R 

X  and  k 
E 
v 


material  cohesion 
material  angle  of  friction 
cap  eccentricity 

slopes  of  the  loading  (compression)  and  unloading  (elastic)  lines 
Young’s  modulus 
Poisson’s  ratio. 


A  good  overview  of  the  development  of  plasticity  theory  and  constitutive 
modeling  of  soil  is  given  in  Scott  (1985).  Schofield  extended  plasticity  theory  to 
the  critical  state  concept,  defining  either  contractile  or  dilatant  deformation  of 
porous  material  as  a  function  of  its  specific  volume  or  void  ratio  (Schofield  and 
Wroth  1968).  In  critical  state  theory,  this  rule  is  developed  around  the  concept  of 
a  “critical  state,”  where  the  plastic  shearing  deformation  occurs  at  a  constant  vol¬ 
ume.  Perhaps  the  most  famous  critical  state  model,  the  Cam-clay  model,  was 
developed  based  on  the  behavior  of  clays.  The  concepts  are  equally  applicable  to 
defining  the  shearing  and  volumetric  behavior  of  granular  materials  such  as 
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granular  soils  or  snow  (Wood  1990).  Although  the  concepts  are  applicable  for 
both  cohesive  and  granular  materials,  the  behavior  of  the  granular  materials  has 
not  been  explored  as  thoroughly,  particularly  regarding  the  influence  of  devia- 
toric  stress  on  the  yield  surface,  which  is  less  clearly  defined  in  soils  but  may 
take  on  a  much  different  shape  than  the  yield  surface  of  metals.  Wood  (1990) 
considers  this  a  difference  in  detail  rather  than  a  difference  in  concept. 

The  modified  CDP  model  also  has  the  features  of  a  critical  state  model  (i.e. 
regions  of  constant  volume  shear  deformation  and  compactive-dilatant  flow). 

The  flow  is  non-associative  (i.e.  the  flow  potential  is  not  associated  with  the  yield 
surface),  except  on  the  cap  surface  of  the  CDP  model. 

The  features  of  the  model  are  best  illustrated  in  terms  of  stress  invariants  for 
effective  stress  and  deviatoric  stress,  p  and  q.  The  first  invariant,  effective  pres¬ 
sure  stress,  p,  is  defined  as 


p  =  --,race(c)=--I,  (]) 

where /j  =  trace  (a)  and  g  is  the  stress  tensor  (the  overscore  denotes  matrix). 

For  a  triaxial  test,  Gi  =  G3  and  eq  1  can  be  reduced  to 


P  =  - 


f  Gj  +  2g3 


(2) 


The  invariant  used  by  ABAQUS,  q,  is  the  deviatoric  stress  defined  as 


q  =  (3) 

where  denotes  the  scalar  product  and 

S  =  g  +  pi  ^ 

where  /  is  the  identity  matrix.  For  a  triaxial  stress  state,  eq  4  is  reduced  to  the 
expression  for  deviatoric  stress: 


9  =  Gt  -  g3 


(5) 


where  g  is  the  total  major  stress  and  g3  is  the  radial  stress  or  confining  pressure. 
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Figure  2.  Modified  Drucker-Prager  yield  sur¬ 
face  in  deviatoric  space.  (After  HKS  2002.) 

For  the  modified  CDP  model  used  in  this  study,  the  yield  surface  is  a  modi¬ 
fied  von  Mises  yield  surface  (i.e.,  the  material  constant  K=  1 .0),  which  is  circular 
in  the  deviatoric  plane  (Fig.  2).  In  the  p-q  plane  the  yield  surface  has  two  major 
segments:  the  Drucker-Prager  portion  of  the  curve  (analogous  to  the  Mohr- 
Coulomb  line)  defines  shear  deformation,  and  the  cap  portion  of  the  surface 
defines  the  intersection  with  the  pressure  axis  and  compressive  deformation, 
shown  in  Figure  3.  The  following  equations  define  the  yield  criteria  in  each  sec¬ 
tion  of  the  yield  surface. 

For  the  Drucker-Prager  shear  or  distortional  failure 


Fs  -  q  -  p  tan  P  -  d  -  0 . 

where  d  is  material  cohesion  and  P  is  material  angle  of  friction. 
For  the  cap  region  of  compactive-dilatant  failure 


^=J(P~PJ2  + 


Rq 


(1  +  a-a/cosP) 


-R[d  +  pa  tanP)  =  0 


(6) 


(7) 


where  a  is  a  transition  parameter,  ranging  typically  from  0.0  to  0.05,  that 
smoothes  the  transition  between  the  shear  failure  and  the  cap  failure;  R  is  a 
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Figure  3.  Modified  Capped  Drucker-Prager  yield  surface  in  the  p-t  plane. 
(After  HKS  2002.) 

material  parameter  controlling  the  cap  eccentricity;  and  p;i  defines  the  evolution 
of  the  cap  along  the  pressure  axes  (via  £>t,-£voi  hardening  law)  according  to 


R  = 


Pb~Pa 


(8) 


Pb  -  Rd 

Pa  ~  (l  +  i?tanP)  ‘ 


(9) 


The  transition  between  the  shear  and  the  cap  failure  is: 


Ft=J(p-pJ  + 


1  — 


a 


cosp 


(d  +  pa  tanp) 


-a(d  +  p&  tanp)  =  0  (10) 


The  pressure-volume  relationships  define  both  hardening  and  softening  through 
volume  changes  based  on  how  the  cap  portion  of  the  yield  surfaces  expands  and 
contracts.  The  hardening  of  the  cap  can  be  defined  as  an  exponential  relationship 
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Figure  4.  Cap  contraction  or  expansion  reflecting  the  softening  (left)  or 
hardening  (right)  of  the  material.  (After  HKS  2002.) 

in  the  pressure-volume  space  as  shown  in  Figure  4,  or  it  is  entered  in  a  piecewise 
fashion  as  a  table  of  and  pairs.  The  piecewise  approach  is  recommended 
by  HKS  (2002)  for  a  better  fit  to  the  data  and  better  model  performance. 

The  main  elements  of  the  modified  CDP  model  are  illustrated  in  Figures  3 
and  4. 

The  plastic  flow  is  defined  by  an  elliptical  potential  surface  as  shown  in  Fig¬ 
ure  5.  Flow  is  associative  (normal  to  the  surface)  in  the  cap  region;  therefore,  the 
equation  for  the  flow  surface  is  identical  to  the  equation  for  the  cap  yield  surface. 
In  the  transition  and  shear  region,  the  flow  is  non-associative  (flow  potential  is 
independent  of  the  failure  surface)  and  is  defined  by  the  following  equation  for 
the  flow  potential: 


Gs=j[(^“-P)tanPl  + 


1  +  a-a/cosp 


(11) 


8 


ERDC/CRREL  TR-05-3 


The  elastic  strains  are  estimated  with  material  parameters  X  and  k,  slopes  of  the 
loading  (compression)  and  unloading  (elastic)  portion  of  the  curve,  respectively, 
E  =  Young’s  modulus,  and  v  =  Poisson’s  ratio. 

Soil  Material 

Triaxial  compression  tests  were  performed  on  Lebanon  Sand.  This  a  highly 
frost-susceptible  material  according  to  laboratory  frost  heave  tests  documented  in 
Shoop  (1988),  is  classified  as  silty  sand  (SM)  using  the  Unified  Soil  Classifica¬ 
tion  System  (USCS),  and  exhibits  thaw  weakening  in  pavement  structures.  This 
soil  would  represent  a  spring  off-road,  trail,  or  subgrade  material  or  a  poor  unsur¬ 
faced  road  material  in  the  northeastern  United  States.  The  soil  was  used  in  full- 
scale  tests  of  paved  and  unpaved  test  sections  in  the  CRREL  Frost  Effects 
Research  Facility  (FERF).  The  material  grain  size  distribution  is  shown  in  Figure 
6.  The  soil  is  non-plastic  (PI  =  0)  and  has  a  maximum  density  of  1826  kg/m3  (114 
pcf)  at  an  optimum  moisture  content  of  13%  (ASTM  D  1557  Modified  Proctor 
Test,  Figure  7).  This  soil  is  characterized  for  moisture  migration  during  freeze- 
thaw  in  Shoop  and  Bigl  (1997),  for  frost  heave  in  Shoop  and  Flenry  (1991),  for 
in-situ  shear  strength  in  Shoop  (1993),  and  for  triaxial  testing  on  frozen  cores  in 
Shoop  (1988).  The  pertinent  background  data  for  Lebanon  Sand  is  included  in 
Appendix  A.  Soil  strength  results  from  California  Bearing  Ratio  (CBR)  tests  at 
various  moisture  contents  are  shown  in  Figure  8  (Shoop  1988). 
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Figure  6.  Lebanon  Sand  grain  size  analysis,  with  D10  =  0.027,  Cu  =  4.579. 
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Figure  7.  Compaction  curve,  with  maximum  density  =  1826  kg/m 
(114  lb/ft3)  and  optimum  moisture  content  =  13%. 


Figure  8.  Soil  strength  from  California  Bearing 
Ratio  (CBR)  tests  at  various  moisture  contents. 
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Sample  Preparation  and  Triaxial  Test  Procedures 

The  tests  performed  on  this  soil  included  saturated  triaxial  tests  and  unsatu¬ 
rated  triaxial  tests  at  specific  moisture  and  density  conditions  experienced  during 
the  freeze-thaw  testing  in  FERF;  the  average  moisture  content  was  13%  (wet  of 
optimum)  and  the  average  dry  density  was  1760  kg/m3  (110  pcf).  Standard  triax¬ 
ial  test  specimens  (2.8  in.  in  diameter  and  5.5  in.  high)  were  fabricated.  The  ini¬ 
tial  and  test  conditions  (moisture  and  density)  for  each  of  the  samples  are  given 
in  Appendix  B. 

All  specimens  were  molded  using  twelve  equal  lifts  of  soil.  Sufficient  water 
was  mixed  with  the  air-dried  soil  to  increase  the  water  content  of  each  lift.  After 
each  lift  was  placed  in  the  mold,  a  handheld  rammer  was  used  to  compact  the  soil 
to  a  specified  thickness.  The  mold  was  then  placed  on  the  base  of  the  triaxial 
chamber  and  placed  under  a  vacuum  of  about  35  kPa  (5  psi).  After  the  vacuum 
had  stabilized  and  it  was  determined  that  the  membrane  was  not  leaking,  the 
mold  was  removed  from  the  specimen,  and  the  dimensions  of  the  membrane- 
encapsulated  specimen  were  recorded.  The  triaxial  chamber  was  then  assembled. 

The  vacuum  applied  to  the  specimen  through  the  pore  pressure  system  was 
then  gradually  reduced  as  the  chamber  pressure  was  simultaneously  increased  to 
the  desired  consolidation  pressure.  The  sample  was  then  back-pressure  saturated, 
maintaining  the  test  effective  stress  until  a  minimum  Skemptom’s  pore  pressure 
parameter  (B)  value  of  0.95  was  obtained  (a  reasonable  value  for  sands).  This 
usually  took  about  a  day  and  resulted  in  a  fully  saturated  sample  (saturation  at  or 
near  100%). 

Once  the  //-value  had  been  achieved,  the  closed-loop  hydraulic  system  was 
activated,  and  the  (axial)  load  cell  was  attached  to  the  specimen’s  top  platen.  Ini¬ 
tial  instrumentation  readings,  i.e.,  chamber  pressure,  pore  pressure,  axial  load, 
and  axial  deformation  using  the  “closed-loop”  LVDT  (linear  variable  differential 
transformer)  and  the  “specimen”  LVDT,  were  recorded.  A  deformation  (or  dis¬ 
placement)  rate-controlled  shear  test  was  conducted.  The  rate  of  strain  was  about 
1%  per  minute.  Following  the  completion  of  each  test,  the  specimen  was 
removed  from  the  triaxial  chamber  and  dried  in  an  oven  at  1 1 0°C  to  obtain  the 
post-test  water  content  and  dry  weight  of  soil. 

Two  types  of  saturated  triaxial  tests  were  conducted:  1)  back-pressure  satu¬ 
rated,  consolidated,  undrained  triaxial  shear  tests  with  pore  pressure  measure¬ 
ments  (R ),  also  known  as  CU  (Confined  Undrained  back-pressure  saturated), 
and  2)  back-pressure  saturated,  consolidated,  drained  triaxial  shear  tests  ( S)  with¬ 
out  pore  pressure  measurements.  In  general  the  procedures  outlined  in  EM 
(Engineer  Manual)  1 110-2-1906  (U.S.  Aimy  COE  1970)  were  followed  for  both 
tests. 
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3  TEST  RESULTS  AND  MATERIAL  PROPERTY 
DETERMINATION 


Saturated  Triaxial  Tests:  Drained  Shear  Result  for  Lebanon  Sand 

Drained  shear  tests  were  performed  using  three  nominal  confining  pressures: 
14,  28,  and  40  kPa  (2,  4,  and  6  psi)  and  two  repetitions  of  each  condition.  The 
low  confining  pressures  were  chosen  to  simulate  the  near-surface  conditions  (for 
a  thawing  off-road  or  trail  surface).  Each  test  was  performed  on  a  7-cm-  (2.8-in.-) 
diameter  by  14-cm-  (5.5-in.-)  high  sample,  which  was  back-pressure  saturated 
and  consolidated.  The  samples  were  prepared  to  initial  conditions  of  9%  water 
content,  approximately  1730  kg/m3  (108  pcf)  dry  density  (compaction  ratio  of 
93%),  and  0.53  void  ratio.  The  water  content  of  the  specimens  after  saturation 
ranged  from  17.5  to  17.8%.  Five  of  these  tests  are  used  in  the  following  analyses 
and  summarized  in  Table  1.  Test  results  are  shown  in  the  two  graphs  in  Figure  9 
as  deviatoric  stress,  q,  vs.  axial  strain,  sax,  and  volumetric  strain  vs.  axial  strain. 

The  shear  failure  envelope  is  obtained  by  plotting  the  data  in  p-q  space  and 
joining  the  yield  points  of  the  stress  paths  of  each  of  the  tests,  as  shown  in  Figure 
10. 

An  alternate  method  of  calculating  the  Drucker-Prager  cohesion,  d  and  //,  is 
to  use  the  following  equations  based  on  Mohr-Coulomb  parameters: 


tan  P  =  1.73  sincp 

(12) 

d  =  cX.l'i  cost)). 

(13) 

Table  1.  Summary  of  Lebanon  Sand  drained  shear  data. 


Test 

no. 

Average 
a3  (kPa) 

Peak 

deviatoric 
stress,  q  (kPa) 

Corresponding 

strain 

Max.  mean 
stress,  p 
(kPa) 

Initial 

void 

ratio 

Max. 

volumetric 

strain 

2DS1 

14 

42 

0.18 

28 

0.53 

0.029 

2DS2 

15 

41 

0.24 

29 

0.54 

0.028 

4DS2 

28 

76 

0.16 

52 

0.53 

0.032 

6DS1 

43 

133 

0.19 

88 

0.54 

0.038 

6DS2 

39 

107 

0.20 

76 

0.54 

0.039 
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Figure  9.  Drained  triaxial  tests. 
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Figure  10.  Shear  failure  line  for  Lebanon 
Sand  drained  shear  tests,  d=  0,  J3=  55.6°. 

The  Mohr-Coulomb  analysis  of  these  tests  is  shown  in  Figure  1 1.  A  failure 
line  drawn  tangent  to  the  failure  circles,  defined  by  eq  14,  results  in  a  cohesion 
value  of  c  =  0  kPa  and  a  friction  angle  of  c|)  =  36°. 

T  =  c  +  otanij).  (14) 

Unfortunately  this  test  series  did  not  include  unloading  and  reloading  the 
sample,  so  the  elastic  properties  were  obtained  by  taking  the  initial  slope  of  the 
stress-strain  curve.  Young’s  modulus,  E,  for  the  drained  shear  tests  is  found  to  be 
17253  kPa.  A  Poisson’s  ratio,  v,  of  0.32  is  found  by  plotting  radial  strain  with 
axial  strain.  It  should  be  noted  that  the  Young’s  modulus  from  the  initial  loading 
portion  may  be  off  by  a  factor  of  two  when  compared  with  that  obtained  from 
unloading-reloading  data  (although  repeatability  in  resilient  modulus  testing  is  of 
the  same  order  of  magnitude).  This  may  not  be  critical  in  our  case  as  we  antici¬ 
pate  the  elastic  strains  to  be  insignificant  compared  to  the  plastic  strains. 
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Figure  11.  Mohr  circle  analysis  for  drained  triaxial  tests,  c  =  0  kPa,  (j>  =  36°. 


Saturated  Triaxial  Tests:  Undrained  Shear  Result  for  Lebanon  Sand 

Undrained  shear  tests,  with  pore  pressure  measurements,  for  Lebanon  Sand 
material  were  also  performed  with  three  confining  pressures:  14,  28,  and  40  kPa 
(2,  4,  and  6  psi)  and  three  test  repetitions.  The  specimens  were  molded  into  a  7- 
cm-  (2.8-in.-)  diameter  by  14-cm-  (5.5-in.-)  high  sample,  then  back-pressure 
saturated  and  consolidated.  The  initial  sample  conditions  were  9%  water  content 
and  approximately  1730  kg/nr'  (108  pcf)  dry  density.  Five  of  these  tests  are  used 
in  the  following  analyses  and  are  summarized  in  Table  2.  A  graph  of  deviatoric 
stress  against  axial  strain  of  the  test  data  is  shown  in  Figure  12.  Pore  pressure 
measurements  are  shown  in  Figure  13. 


Table  2.  Summary  of  Lebanon  Sand  undrained  shear  data. 


Test 

no. 

Average  cr3 
(confining 
pressure)  (kPa) 

Peak  q 
(<*3  ~  cr3) 
(kPa) 

Corresponding 

Initial  void 
ratio,  e 

a'3 

(kPa) 

(kPa) 

P' 

(kPa) 

2US3 

13 

25 

9 

35 

22 

0.52 

4US2 

28 

21 

13 

34 

24 

0.51 

4US3 

26 

20 

10 

30 

20 

0.53 

6US2 

40 

26 

17 

45 

31 

0.53 

6US3 

40 

27 

17 

44 

31 

0.50 
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Figure  12.  Consolidated,  undrained  triaxial  tests. 
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Figure  13.  Undrained  triaxial  tests  pore  pressure  measurements. 
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The  stress  paths  from 
each  test,  in  p'-q'  space, 
are  shown  in  Figure  14. 

As  with  the  drained  data, 
the  critical  state  analysis 
for  undrained  shear  data 
was  derived  from  eq  2 
and  5  using  the  effective 
stresses.  A  compilation  of 
the  test  data  for  the 
saturated  conditions  is 
given  in  Figure  15  and 
shows  the  good  agree¬ 
ment  among  the  tests. 

The  Mohr’s  circle 
analysis  of  these  tests  is 
shown  in  Figure  16.  A 
failure  line  is  drawn  tan¬ 
gent  to  the  remaining  four 
failure  circles  having  a 
cohesion  value  of  c  =  4.8 
kPa  and  a  friction  angle 
of  ((>  =  17°  for  effective 
stress,  and  c=  12.1  kPa  Figure  15.  Test  data  for  all  of  the  saturated 

and  (()  =  0°  for  total  stress.  triaxial  tests  (drained  and  undrained). 
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Figure  16.  Mohr  circle  analysis  for  undrained  triaxial  tests;  c  =  12.1  kPa  and  c|>  =  0°  for  total 
stress  (top),  and  c  =  4.8  kPa  and  cj>  =  17°  for  effective  stress  (bottom). 


Hydrostatic  Consolidation  Test 

Hydrostatic  consolidation  tests  were  performed  to  determine  the  pressure- 
volume  relationship  of  the  material.  The  test  results  are  shown  graphically  in 
Figure  17  and  are  compared  to  a  similar  sand  in  a  compacted  state  in  Figure  18. 
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Volumetric  Strain 


Figure  17.  Hydrostatic  compression  vs.  volumetric  strain. 


Figure  18.  Pressure-volume  relationship  of  Lebanon  Sand  compared  to 
compacted  sand. 
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The  pressure-volume  relationship  can  also  be  input  directly  in  tabular  form. 
The  latter  is  suggested  as  more  accurate  (HKS  2002)  and  was  used  in  this  study 
(Table  3). 


Table  3.  Hardening  model  for  Lebanon  Sand. 


Hydrostatic  stress,  Pb 

MPa  (psi) 

Volumetric  plastic 
strain 

0.0082  (1.2) 

0 

0.0389  (5.6) 

0.009 

0.0760  (11.0) 

0.022 

0.1639  (23.8) 

0.038 

0.3655  (53.0) 

0.054 

0.7201  (104.4) 

0.072 

Unloading 

Volumetric  strain 

0.3569  (51.8) 

0.071 

0.0822  (11.9) 

0.067 

0.0178  (2.6) 

0.064 

Constant  Water  Content  Tests  at  Specific  Thawing  Conditions 

Several  triaxial  tests,  without  pore  pressure  measurements,  were  performed 
on  Lebanon  Sand  at  an  average  dry  density  of  1760  kg/m3  (110  pcf)  and  moisture 
content  of  13%  (68%  of  saturation).  A  brief  summary  of  the  results  are  presented 
here.  Additional  analysis  and  a  more  detailed  report  on  these  tests  are  forthcom¬ 
ing.  The  specimens  were  molded  into  a  7-cm-  (2.8-in.-)  diameter  by  15-cm-  (6- 
in.-)  high  cylinders  in  five  equal  lifts.  The  confining  pressures  used  were  14,  28, 
and  40  kPa  (2,  4,  and  6  psi).  Each  test  was  subjected  to  loading  and  unloading 
cycles  as  shown  in  Figure  19.  Table  4  summarizes  the  test  analysis. 


Table  4.  Summary  of  Lebanon  Sand  undrained 
shear  tests  for  samples  at  13%  moisture  content. 


Test 

no. 

Average  cr3 
(confining 
pressure)  (kPa) 

Peak  q 
(CTi  -  ct3) 
(kPa) 

Max.  mean 
stress,  p 
(kPa) 

2.1 

14 

57 

33 

2.3 

14 

61 

34 

2.4 

14 

86 

43 

2.5 

14 

101 

47 

4.2 

27 

83 

55 

6.1 

41 

159 

94 

6.2 

42 

121 

82 
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Figure  19.  Triaxial,  undrained  shear  data  for  samples  representing  specific 
thawing  conditions  of  an  average  moisture  content  of  13%  and  a  dry  density  of 
1760  kg/m3.  Legend  indicates  test  number  from  Table  4  and  confining  pressure  in  psi. 


The  p-q  yield  analyses 
resulted  in  d  =  16.7  kPa,  P  = 
54.8°  (Fig.  20).  These  tests 
show  some  cohesion  value,  or 
apparent  cohesion,  as  is  often 
the  case  in  unsaturated,  moist 
soil. 

The  Mohr’s  circle  analy¬ 
sis  is  given  in  Figure  21.  The 
resulting  c  and  (p  values  are 
very  similar  to  those  obtained 
from  the  saturated,  drained 
tests. 

A  Young’s  modulus 
value  of  8490  kPa  was 
obtained  by  averaging  all  of 
the  values  taken  from  the 
stress  vs.  strain  graphs  of 
each  laboratory  test.  The 
slope  of  each  loading  curve 
was  measured  and  then 


Figure  20.  Shear  failure  line  for  samples  at 
13%  moisture  content,  d=  16.7  kPa,  p  =  54.8°. 
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Figure  21.  Mohr’s  circle  analysis  for  samples  at  13%  moisture  content;  c  =  9.4  kPa,  (|>  =  33°. 

averaged  with  the  other  loading  curves  for  that  specific  graph.  An  average  was 
found  for  each  confining  pressure  [14,  28,  and  40  kPa  (2,  4,  and  6  psi)]  and  then 
a  final  Young’s  modulus  of  8490  kPa  (1230  psi)  was  calculated.  Graphs  of  the 
test  data  and  a  table  of  the  calculations  are  given  in  Appendix  C. 

A  summary  of  the  shear  failure  line  parameters  obtained  using  the  Mohr- 
Coulomb  analysis  and  eq  12  and  13  are  given  in  Table  5.  In  Table  5,  d  ranges 
from  0  to  13  kPa,  and  P  ranges  from  27  to  45°.  These  are  somewhat  lower  than 
the  values  from  the  p-q  plots  and  were  not  used  for  the  remainder  of  the  study. 


Table  5.  Shear  line 
Coulomb  analysis. 

parameters 

from 

Mohr- 

c 

(kPa) 

d 

(kPa) 

<l> 

n 

P 

n 

Drained,  saturated 

0 

0 

36 

45 

Undrained,  saturated 

4.8 

8 

17 

27 

Undrained,  13%  me 

9.0 

13 

33 

44 

Summary  of  Test  Data 

Table  6  summarizes  the  material  parameters  from  all  of  the  test  data.  In  gen¬ 
eral,  the  constitutive  model  parameters  shown  in  Table  6  reflect  the  specific  test 
conditions.  The  final  model  parameters  chosen  to  represent  the  thawing  soil 
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behavior  were  a  combination  of  the  parameters  obtained  from  the  different  test 
conditions.  The  Young’s  modulus  values  from  the  loading-reloading  curve  were 
about  half  of  those  from  the  initial  loading  of  the  saturated  triaxial  tests.  As  most 
of  the  deformation  is  plastic,  the  value  of  Young’s  modulus  does  not  have  a  huge 
influence  on  the  results,  and  the  value  from  the  unsaturated  samples  was  chosen 
to  represent  the  soil  in  the  final  model  parameters.  The  other  parameter  that  was 
significantly  between  the  saturated  tests  and  the  tests  at  13%  moisture  content  is 
the  cohesion.  Both  of  the  saturated  test  conditions  show  essentially  no  cohesion, 
as  would  be  expected  for  saturated,  sandy  soil.  The  samples  tested  at  13%  mois¬ 
ture  content  show  a  small  amount  of  cohesion,  or  apparent  cohesion,  as  is  typical 
in  unsaturated  soils  because  of  moisture  tension.  A  value  between  these  two 
extremes  was  chosen  to  represent  the  overall  thawing  soil  behavior  in  the  final 
constitutive  model.  The  remaining  parameters  were  within  the  range  of  the  tests 
and  were  refined  to  obtain  the  best  response  in  model  simulations  as  discussed  in 
the  following  section. 

The  p  values  appear  high,  but  they  are  based  on  the  ABAQUS  formulation  of 
the  deviatoric  stress  invariant  q  and  agree  with  values  reported  by  Arnold  (2002). 
Using  a  Mohr-Coulomb  approach,  the  angle  of  internal  friction  cp  for  these  tests 
is  36°  and  17°  for  saturated  drained  and  undrained  tests,  respectively,  and  33°  for 
the  tests  at  13%  moisture  content;  these  values  are  reasonable  for  sands. 


Table  6.  Material  constants  comparison  for  the  Capped  Drucker-Prager  model  on  Lebanon 
Sand. 


Material  parameter 

Saturated 
drained  shear 

Saturated 
undrained  shear 

Shear  at  13% 
moisture  content 

Final  model 
parameters 

E,  Young’s  modulus  (kPa) 

17,250 

17,440 

8,490 

8500 

v,  Poisson’s  ratio 

0.32 

0.5 

— 

0.32 

P,  Drucker-Prager  angle 
of  friction  (degrees) 

55.6 

56.4 

54.8 

55.8 

d,  Drucker-Prager 
material  cohesion 
(kPa) 

0 

0 

16.7 

10.0 

R,  cap  eccentricity 

0.15-0.88 

0.45 

Initial  value  of  volumetric 
plastic  strain 

0.001 

0 

K,  flow  stress  ratio 

1.0 

1.0 
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4  SIMULATION  OF  LABORATORY  TESTING 

The  material  properties  generated  from  laboratory  testing  (Table  6)  were 
used  to  create  a  working  simulation  of  the  biaxial  load  tests  by  manipulating  the 
parameters  of  the  material  model  to  ensure  a  good  fit  for  a  wide  range  of  simu¬ 
lated  test  parameters.  After  a  good  fit  to  the  triaxial  test  data  was  obtained,  the 
material  model  was  used  in  a  3-D  simulation  of  an  unpaved  road.  Some  addi¬ 
tional  adjustments  to  the  material  model  were  needed,  primarily  to  account  for 
apparent  cohesion  in  the  unconfined  road  surface.  Although  the  best- fit  parame¬ 
ters  to  the  saturated  biaxial  test  data  included  very  low  values  of  CDP  cohesion 
(i d  ~  0),  this  resulted  in  immediate  dilatant  failure  of  the  surface  elements  as  the 
wheel  was  applied  to  the  soil.  Arguably  this  surface  layer  is  likely  not  saturated 
and  some  apparent  cohesion  would  be  present  (as  in  the  13%  moisture  content 
test  results),  so  some  cohesion  was  included  in  the  final  constitutive  model 
parameters. 

A  final  set  of  parameters  was  obtained  that  gave  a  reasonable  fit  to  both  lab 
and  field  observations.  These  parameters  are  listed  in  the  last  column  of  Table  6. 
The  final  material  model  parameters  are  compared  to  the  laboratory  data  in  Fig¬ 
ure  22.  The  model  captures  the  stiffness  of  the  higher  confining  pressure  better 
than  that  of  the  low  confining  pressure. 


o.oo 
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o.io 


0.15 

Axial  Strain 


0.20 


0.25 


0.30 


Figure  22.  Comparison  of  final  model  parameters  (Table  6)  with  laboratory 
data  for  14,  28,  and  40  kPa  confining  pressure. 
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5  VALIDATION  AND  APPLICATION 

Finite  element  simulations  of  experiments  performed  on  Lebanon  Sand  were 
used  to  validate  the  material  model  and  assure  that  it  was  robust  under  a  variety 
of  loading  and  simulation  conditions.  Three  finite  element  simulations  were  con¬ 
structed;  details  of  these  finite  element  simulations  are  given  elsewhere.  For 
completeness  of  this  report,  brief  descriptions  of  each  simulation  are  listed  below. 

•  Direct  shear  tests  were  performed  on  thawing  Lebanon  Sand  samples 
taken  from  the  FERF  test  sections.  The  experimental  data  for  the  direct 
shear  tests  are  summarized  in  Shoop  (1993).  Tests  were  completed  at  six 
values  of  normal  loading.  The  mesh  for  the  simulations  is  shown  in  Fig¬ 
ure  23a.  The  simulations  modeled  the  data  very  well  except  for  very 
small  normal  loads,  as  shown  in  Figure  24  (Flaehnel  and  Shoop,  in 
prep.). 

•  Various  finite  element  models  were  generated  to  simulate  vehicle  travel 
on  an  unpaved  road  or  trail  or  an  off-road  condition.  These  models  were 
implemented  in  a  three-dimensional  dynamic  analysis  with  a  wheel  roll¬ 
ing  over  the  soil.  Simulations  of  the  unpaved  road  were  used  to  study  the 
effects  of  vehicle  speed,  load,  suspension  system,  wheel  torque,  wheel 
slip,  etc.,  on  the  deterioration  of  unsurfaced  roads  or  trails  in  terms  of 
rutting  and  washboard  formation  (Shoop  et  al.  2002).  One  of  the  meshes 
for  this  study  is  shown  in  Figure  23b.  Qualitatively  the  results  from  the 
finite  element  model  agree  with  observations  of  rutting  and  approximate 
resisting  forces  on  the  wheel.  Both  the  simulation  and  field  observations 
indicate  moderate  plowing  and  flowing  of  the  soil  to  the  sides  and  in 
front  of  the  wheel  and  the  wheel  sinking  through  the  full  depth  of  the 
thawed  soil  layer  such  that  the  wheel  is  nearly  rolling  on  the  top  of  the 
frozen  layer. 

•  The  final  study  was  for  a  layered  pavement  system  (Fig.  23c).  This 
model  will  be  used  with  varying  pavement  structures  and  frost-thaw 
depths  to  study  the  impact  of  different  vehicle  loading  and  frost-thaw 
layering  on  the  pavement  deformation.  This  type  of  model  will  also  be 
used  to  study  airfield  structures  for  unpaved  and  marginal  airfields  dur¬ 
ing  thaw  or  for  paved  airfields  with  thawing  subgrades. 
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a.  Direct  shear  test. 


c.  Layered  pavement  system. 

Figure  23.  Three-dimensional  finite  element  simulations,  all  incorporating  the  Lebanon 
Sand  thawing  soil  material  model. 
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Figure  24.  Comparison  of  finite  element  simulations  of  the  direct  shear  test  with  laboratory 
data. 
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6  THAWING  SOIL  MODEL  COMPARED  TO  OTHER 
GRANULAR  MATERIALS 

The  CDP  model  for  Lebanon  Sand  is  compared  to  some  other  granular  mate¬ 
rials  in  Table  7.  The  McCormick  Ranch  Sand  (HKS  1996,  DiMaggio  and  Sandler 
1971,  1976)  and  the  Lebanon  Sand  are  fairly  similar  in  grain  size,  while  the  A-4 
soil  (Burton  et  al.  2003)  is  a  silt  (ML).  The  grain  size  distributions  are  plotted  in 
Figure  25. 

Figure  26  compares  the  hardening  curves  for  the  materials.  Although  the 
McCormick  Ranch  Sand  has  a  grain  size  distribution  similar  to  that  of  Lebanon 
Sand,  the  McCormick  Sand  was  tested  in  a  much  more  compacted  state  and  at 
higher  pressures  (for  a  much  different  application).  Also,  the  A-4  soil  samples 
were  prepared  at  optimum  moisture  and  density,  not  the  loose  conditions  of  the 
thawing  Lebanon  Sand.  Therefore,  both  other  soils  required  higher  pressures  than 
Lebanon  Sand  for  the  same  volumetric  strain. 

For  contrast,  a  fresh  snow  material  model  from  Shoop  (200 1 )  was  also 
included.  The  compaction  of  snow  is  an  order  of  magnitude  greater  than  for  soils, 
as  can  be  seen  by  the  scale  for  volumetric  strain  of  0  to  1 .6  for  snow  and  0  to  0. 1 
for  the  soils  (Fig.  26). 


Table  7.  Material  constants  for  Capped  Drucker-Prager  model  for  various  materials. 


Material  parameter 

A-4*  saturated 
drained  shear 

McCormick  Ranch 
Sand*  saturated 
drained  shear 

Snow** 

(SNOW2) 

Lebanon 
Sand  final 

£,  Young’s  modulus,  kPa 

16,000 

689,475 

1,379 

8,500 

v,  Poisson’s  ratio 

0.33 

0.25 

0.3 

0.32 

P,  Drucker-Prager  angle  of 
friction 

53° 

14.6 

22.5 

55.8 

d,  Drucker-Prager  material 
cohesion,  kPa 

16.6 

1190 

15 

10 

R,  cap  eccentricity 

0.45 

0.1 

1.1e-4 

0.45 

s  ,pl  L. ,  initial  value  of 

0.001 

0.001 

0.001 

0 

volumetric  plastic  strain 

K,  flow  stress  ratio 

1.0 

1.0 

1.0 

1.0 

*  Burton  et  al.  (2003) 

f  HKS  (1996),  DiMaggio  and  Sandler  (1971,  1976) 
**  Shoop  et  al.  (1999),  Shoop  (2001) 
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Figure  25.  Lebanon  Sand  grain  size  analysis,  D10  =  0.027,  Cu  =  4.579, 
compared  with  the  McCormick  Sand  and  A-4  silt. 
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Figure  26.  Pressure-volume  relationship  of  Lebanon  Sand 
compared  with  some  other  granular  materials. 
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Figure  27  compares  the  shear  failure  lines  for  the  thawing  Lebanon  Sand  at 
three  test  conditions,  the  A4  silt,  and  snow.  The  shear  lines  are  surprisingly  dif¬ 
ferent  for  the  McCormick  Sand  (not  shown  in  Fig.  27),  which  plots  far  above  the 
other  materials  because  of  its  large  cohesion  value.  The  slope,  p,  values  cover  a 
wide  range  also  (14.6°  to  55.8°),  with  the  snow  falling  in  the  middle  of  the  range 
at  22.4°. 


Figure  27.  Drucker-Prager  shear  failure  comparison  of  Lebanon  Sand 
tested  at  different  conditions,  A-4  silt  (saturated,  drained),  and  snow. 
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7  SUMMARY  AND  CONCLUSIONS 

Spring  thaw  creates  unique  conditions  for  soil.  The  thawing  soil  is  often 
abnormally  loose  and  wet  as  a  result  of  the  expansion  from  freezing  and  the  sub¬ 
sequent  excess  meltwater  being  trapped  in  the  thawing  later  by  impermeable  fro¬ 
zen  layers  below.  The  thaw  layer  is  very  weak,  and  it  deforms  plastically  in  both 
compaction  and  shear. 

A  frost-susceptible  fine  sand,  Lebanon  Sand,  was  subjected  to  a  range  of  tri- 
axial  testing  using  density,  moisture,  and  loading  conditions  similar  to  those 
experienced  during  the  freeze-thaw  testing  in  CRREL’s  Frost  Effect  Research 
Facility.  A  Capped  Drucker-Prager  (CDP)  plasticity  model  was  used  to  simulate 
the  thawed  soil  using  the  commercial  finite  element  code  ABAQUS.  Simulations 
of  the  triaxial  tests  were  compared  to  the  laboratory  data.  The  material  model  was 
then  implemented  in  three-dimensional  finite  element  simulations  to  check  model 
robustness  and  to  validate  it  using  direct  shear  test  data,  experiments  on  vehicle 
mobility  on  thawing  soil,  and  a  layered  pavement  system  subjected  to  freeze- 
thaw. 

The  following  summarize  significant  results  of  this  study: 

•  A  constitutive  model  of  thawing  soil  was  generated  based  on  triaxial  and 
hydrostatic  test  data. 

•  The  model  was  implemented  in  several  3D  finite  models  using  the 
commercial  finite  element  code  ABAQUS. 

•  A  small  apparent  cohesion  is  needed  in  the  material  model  to  obtain 
realistic  behavior  in  the  unconfined  3D  model  structural  model  of  the 
wheel  rolling  on  thawing  ground. 

•  The  3D  simulations  attested  to  the  robustness  of  the  material  model 
through  the  attainment  of  realistic  model  results  under  a  variety  of  con¬ 
ditions. 
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APPENDIX  A.  GEOTECHNICAL  CHARACTERIZATION  DATA 
SUMMARY  FOR  LEBANON  SAND 

A  controlled  full-scale  freeze-thaw  test  using  Lebanon  Sand  was  conducted 
in  the  test  basin  of  the  Frost  Effects  Research  Facility  (FERF).  The  test  area  was 
36.6  x  13.1  m  and  consisted  of  a  1.1-m  layer  of  Lebanon  Sand  and  a  base  layer 
of  fill  separated  by  a  geotextile.  The  soil  was  compacted  to  the  desired  density 
and  at  a  specified  moisture  content.  Temperature  and  soil  moisture  sensors  were 
installed.  Several  freeze-thaw  cycles  were  created,  and  the  moisture  content  in 
the  soil  was  varied  by  changing  the  water  table  depth  or  by  adding  surface  water 
to  simulate  rain  or  snowmelt.  The  soil  is  frozen  from  the  surface  downward  using 
freezing  panels  and  then  thawed  with  an  elevated,  constant  air  temperature.  When 
the  frost  reached  the  required  depth,  frost  heave,  temperature,  and  soil  moisture 
were  measured.  Frozen  cores  were  drilled  during  different  freeze  cycles  to 
determine  ice  contents.  Various  strength  tests  were  performed  in  the  laboratory 
from  frozen  core  specimens  (Shoop  1990,  1993).  When  the  soil  thawed  to  the 
desired  depth,  moisture  content,  density,  and  in-situ  strength  were  measured.  In 
addition,  mobility  testing  was  conducted  using  the  CRREL  Instrument  Vehicle 
(CIV)  to  determine  traction  and  motion  resistance  under  thawing  conditions  on 
half  of  the  test  area.  Physical  and  hydraulic  properties  under  frozen  and  unfrozen 
conditions  were  compiled  (Shoop  and  Bigl  1997).  The  soil  moisture  characteris¬ 
tic  was  determined  in  the  laboratory.  Frost  heave  and  water  migration  during 
freezing  and  thawing  were  characterized  as  well  (Shoop  and  Henry  1991). 

This  appendix  conatins  the  following  sections: 

Al.  Strength  of  Frozen  Samples 

A2.  Physical  and  Hydraulic  Properties  of  Lebanon  Sand 

A3.  Water  Migration  and  Frost  Heave  Characteristics 

A4.  Shear  Measurements 

A5.  Chemical  Analysis. 
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A1.  Strength  of  Frozen  Samples 

Strength  tests  were  conducted  from  frozen  cores  taken  from  the  FERF 
(Shoop  f  988).  Tests  S-f ,  S-5,  and  S-6  were  conducted  at  a  frozen  state  with  zero 
confining  pressure,  while  tests  S-2,  S-3  and  S-4  were  conducted  at  room  tem¬ 
perature  with  a  5-psi  confining  pressure.  The  strain  rate  for  each  test  was  main¬ 
tained  at  approximately  1 0%  strain  per  second,  and  radial  strain  was  measured 
for  the  frozen  samples  only.  Results  for  the  stress-strain  and  Poisson’s  stress 
graphs  are  shown  in  Figures  A.  1 . 1  and  A.  1 .2,  respectively. 


Table  A1.  Summary  for  strength  tests  of  frozen  samples. 


Test  samples 

S-1 

S-2 

S-3 

S-4 

S-5 

S-6 

Peak  force  (lb) 

3516 

54 

71 

32 

4691 

2979 

Peak  stress  (psi) 

1117.9 

16.8 

17.8 

8.1 

1446.4 

933.6 

Strain  at  peak 
stress  (%) 

1.58 

4.06 

23.19 

21.55 

2.74 

1.60 

Displacement  at 
peak  stress  (in.) 

0.079 

0.194 

0.983 

0.989 

0.127 

0.080 

Secant  modulus 
(kpsi) 

70620 

413.15 

76.71 

37.56 

52780 

58300 

Temperature  (°C) 

-4 

-4.4 

-4.5 

Average  strain  rate 
(%/s) 

10.82 

10.71 

10.52 

10.09 

11.34 

10.29 

Poison’s  ratio  at 
half  peak  stress 

0.062 

0.069 

0.056 

Poison’s  ratio  at 
peak  stress 

0.202 

0.214 

0.282 

Strain  rate  at  peak 
stress  (%/s) 

10.02 

10.36 

11.21 

10.84 

10.71 

9.31 

Initial  density  (lb/ft3) 

117.87 

116.59 

114.40 

107.42 

117.72 

113.81 

Time  to  peak  (ms) 

158.0 

392.0 

2068.0 

1988.0 

256.0 

172.0 

Length  of  test  (ms) 

2002.0 

2024.0 

2120.0 

2036.0 

2016.0 

2016.0 

Sample  weight  (g) 

479.19 

461.48 

451.95 

452.17 

458.32 

472.78 

Sample  height  (in.) 

5.002 

4.875 

4.750 

5.102 

4.700 

5.040 

Sample  diameter 
(in.) 

1.986 

1.985 

2.009 

2.001 

2.005 

2.000 

Water  content  (%) 

10.2 

10.1 

9.4 

7.4 

11.6 

9.7 

Test  temperature 
(°C) 

-6.0 

22.0 

22.0 

22.0 

-6.0 

-6.0 

Confining  pressure 
(kPa) 

0 

5 

5 

5 

0 

0 
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Figure  A1.  Stress-strain  plots  from  the  frozen  cores. 
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A2.  Hydraulic  Properties  for  Lebanon  Sand 

The  soil  moisture  characteristics  and  hydraulic  properties  from  the  tests  in 
the  FERF  are  summarized  below  (Shoop  and  Bigl  1997). 


Table  A2.  Physical  and  hydraulic  properties  of  Lebanon  Sand  in  the  FERF  test  sections. 


Properties 

Layer  1 
(0-15  cm) 

Layer  2 
(15-46  cm) 

Layer  3 
(46-107  cm) 

Soil  density  (g/cm3) 

1.551* 

1.678 

1.697 

Soil  porosity  (cm3/cm3) 

0.419 

0.419 

0.419 

Soil  water  characteristics:  Gardner’s  multiplier  for 

moisture,  A„ 

1.962  x  10~5 

1.962  x  10~5 

1.962  x  10~5 

Gardner’s  multiplier  for 
moisture,  a 

1.975 

1.975 

1.975 

Min.  unfrozen  water  cont.  (cm3/cm3) 

0.03602 

0.03602 

0.03602 

Saturated  hydraulic  conductivity  (cm/hr) 

1.6 

1.6 

1.6 

Permeability  Gardner’s  multiplier  for 

characteristics:  hydraulic  cond.,  AK 

1.590  x  10'9 

1.590  x  10'9 

1.590  x  10~9 

Gardner’s  multiplier  for 
hydraulic  cond.,  p 

4.623 

4.623 

4.623 

‘1.551  g/cm3for  dry  case,  1.620  g/cm3forwet  case. 


Figure  A3.  Hydraulic  conductivity  and  moisture  retention  for 
Lebanon  Sand. 
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Table  A3.  Warming  data  from  an  unfrozen  water 
content  test.  Four  samples  tested  with  thawed 
gravimetric  water  contents  are  shown. 


Temp  Unfrozen  water  content  (%), 

(°C)  gravimetric 


>0 

5.52 

10.87 

15.74 

20.65 

-3.00 

0.32 

0.2 

0.18 

-1.93 

0.63 

0.41 

0.39 

0.37 

-1.46 

0.95 

0.62 

0.59 

0.55 

-0.96 

1.26 

0.82 

0.79 

0.74 

-0.78 

1.26 

1.03 

0.99 

0.74 

-0.59 

1.58 

1.03 

0.99 

0.92 

-0.48 

1.89 

1.24 

1.18 

1.29 

-0.37 

1.89 

1.65 

1.38 

1.29 

-0.29 

2.2 

1.65 

1.58 

1.29 

-0.17 

2.52 

2.06 

1.78 

1.67 

-0.08 

4.09 

2.89 

2.57 

2.04 

-0.04 

4.72 

3.51 

2.96 

2.59 

-0.01 

7.64 

6.32 

5.93 

Table  A4.  Coefficients  for 


modeling  unfrozen  water  con¬ 
tent. 

Sample 

water 

content 

Coefficients 

(%) 

a 

P 

5.52 

1.0045 

-0.5629 

10.87 

0.7949 

-0.5069 

15.74 

0.6685 

-0.5297 

20.65 

0.6075 

-0.5235 

Figure  A4.  Warming  data  from  unfrozen  water 
content  test  on  four  samples  tested  with 
thawed  gravimetric  water  contents. 
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A3.  Water  Migration  and  Frost  Heave  Characteristics 

Results  from  the  laboratory  frost  heave  and  moisture  tension  tests  are  shown 
in  Figures  A5  and  A6,  respectively  (Shoop  and  Henry  1991). 


Figure  A5.  Laboratory  frost  heave  test  data. 
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Figure  A8.  Lebanon  Sand  (red  star)  position  on  frost-susceptibility  chart. 
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A4.  Shear  Measurements 

The  standard  laboratory  tests,  including  triaxial  compression  and  direct  shear, 
were  performed  along  with  in-situ  tests  using  a  shear  annulus  and  CIV  (Shoop 
1990,  1993).  Table  A5  shows  the  summary  from  the  triaxial  compression  test. 

The  specimens  used  for  the  triaxial  compression  test  were  all  machined  from  the 
various  frozen  cores  taken  from  the  FERF  test  section.  The  frozen  samples  were 
prepared  for  testing,  then  allowed  to  thaw  overnight,  and  tests  were  performed  at 
room  temperature.  Table  A6  summarizes  the  direct  shear  test  that  was  performed 
from  some  samples  that  were  taken  from  the  FERF  test  section  as  it  thawed  and 
from  other  samples  that  were  prepared  from  the  frozen  cores  similar  to  the  triax¬ 
ial  test  sample  preparation.  The  shear  annulus  test  was  conducted  in  the  FERF 
test  section  immediately  prior  to  the  mobility  tests  with  the  CIV  and  is  listed  in 
Table  A7.  From  the  mobility  tests  using  the  CIV,  the  maximum  shear  and  normal 
stresses  were  developed  from  the  tire-soil  interface.  The  maximum  shear  stress 
applied  to  the  soil  was  derived  from  the  peak  gross  tractive  force  divided  by  the 
contact  area  of  the  tire.  The  normal  stress  was  measured  by  dividing  the  normal 
load  on  the  tire  by  the  tire  contact  area.  The  maximum  shear  and  normal  stresses 
express  in  c'  and  ())'  for  various  soil  conditions  are  listed  in  Table  A8.  The  loading 
comparison  for  each  test  method  are  listed  in  Table  A9  and  graphed  in  Figure 
A9.  The  effect  of  moisture  content  on  (j)'  measured  using  the  direct  shear,  shear 
annulus,  and  CIV  methods  showed  that  (j)'  increased  until  the  soil  liquid  limit  was 
met  and  then  rapidly  dropped  (Figure  A10).  The  calculated  ())'  from  the  CIV  test 
with  thaw  depth  is  also  shown  in  Figure  All. 


Table  A5.  Triaxial  test  data  (thawed  cores). 


Confining 

Failure 

Water 

Dry 

Total 

Sample 
no . 

pressure 

(kPa) 

stress 

(kPa) 

content 

(%) 

density 

(g/cm3) 

density 

(g/cm3) 

s-9 

0 

70.97 

16.0 

1.592 

1.846 

s-6 

0 

163.98 

16.5 

1.672 

1.948 

s-1 

0 

175.01 

18.7 

1.654 

1.963 

s-3 

0 

207.39 

18.8 

1.659 

1.971 

s-4 

34.4 

55.81 

11.6 

1.539 

1.718 

s-3 

34.4 

122.64 

7.4 

1.704 

1.830 

s-2 

34.4 

115.75 

9.7 

1.701 

1.866 

s-5 

68.9 

228.75 

21.7 

1.558 

1.897 

s-12 

68.9 

289.38 

18.4 

1.630 

1.931 

s-7 

137.8 

319.70 

15.1 

1.574 

1.811 

s-11 

137.8 

266.64 

12.5 

1.629 

1.833 

s-4 

137.8 

320.39 

18.6 

1.650 

1.957 

_ 
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Table  A6.  Direct  shear  test  data. 


Sample 

no. 

Normal 

load 

(N) 

Peak 

load 

(N) 

Residual 

load 

(N) 

Area 
< cm ?) 

Water 

content 

(%) 

Dry 
density 
( glcm 3) 

Total 

density 

(glcm3) 

12 

243.48 

146.78 

93.41 

30.53 

21.5 

1.501 

1.823 

8 

243.48 

179.25 

160.13 

30.61 

12.6 

1.573 

1.772 

tcl2r 

243.48 

224.18 

164.58 

31.41 

12.5 

1.696 

1.909 

tcl2 

243.48 

231.30 

157.90 

31.41 

12.3 

1.696 

1.905 

tell 

243.48 

232.63 

151.23 

31.41 

12.6 

1.696 

1.910 

2 

243.48 

236.19 

179.25 

30.64 

10.6 

1.702 

1.883 

7 

360.29 

293.57 

265.10 

30.51 

12.4 

1.675 

1.882 

tcl2 

360.29 

320.26 

271.33 

31.41 

11.6 

1.696 

1.894 

5 

360.29 

324.70 

240.19 

30.66 

18.1 

1.695 

2.002 

tell 

360.29 

326.04 

233.52 

31.41 

12.5 

1.696 

1.909 

tcl2r 

360.29 

346.94 

263.32 

31.41 

12.3 

1.696 

1.905 

4 

360.29 

373.63 

272.22 

30.35 

16.9 

1.735 

2.028 

6 

556.00 

389.20 

345.61 

30.75 

20.5 

1.470 

1.772 

11 

556.00 

411.44 

350.06 

30.51 

10.6 

1.631 

1.804 

9 

556.00 

422.56 

341.16 

30.53 

19.1 

1.614 

1.922 

tcl2r 

556.00 

428.34 

311.36 

31.41 

13.5 

1.696 

1.925 

tcl2 

556.00 

437.02 

311.36 

31.41 

11.9 

1.696 

1.898 

tell 

556.00 

444.80 

400.32 

31.41 

12.6 

1.696 

1.911 

10 

711.68 

494.84 

422.56 

30.41 

21.6 

1.451 

1.765 

tcl2 

711.68 

533.76 

355.84 

31.41 

15.8 

1.696 

1.964 

13 

711.68 

550.44 

442.58 

30.25 

19.7 

1.627 

1.947 

tcl2r 

711.68 

551.55 

389.20 

31.41 

11.5 

1.696 

1.892 

tcl2r 

711.68 

578.24 

378.08 

31.41 

12.7 

1.696 

1.911 

1 

711.68 

633.84 

456.81 

30.57 

17.0 

1.719 

2.011 

Table  A7.  Shear  annulus  test  data. 


Date 

c 

(kPa) 

4> 

(degrees) 

Thaw 

depth 

(cm) 

Water 

content 

(%) 

Dry 

density 

(g/cm3) 

Total 

density 

(glen?) 

Saturation 

<%) 

Nov  9 

4.44 

33.5 

3.81 

20.0 

1.509 

1.811 

67.9 

Nov  10 

0.88 

36.5 

10.16 

23.4 

1.546 

1.907 

83.9 

Nov  10 

0.30 

37.9 

10.16 

23.4 

1.546 

1.907 

83.9 

Nov  10 

0.92 

36.6 

10.16 

23.4 

1.546 

1.907 

83.9 

Dec  20 

0.90 

38.4 

* 

13.0 

1.699 

1.920 

58.9 

Dec  20 

1.54 

37.9 

* 

13.0 

1.699 

1.920 

58.9 

Dec  27 

2.82 

36.8 

* 

9.0 

1.611 

1.756 

35.6 

Dec  27 

4.13 

34.8 

* 

9.0 

1.611 

1.756 

35.6 

Dec  27 

3.17 

33.5 

* 

9.0 

1.611 

1.756 

35.6 

Dec  28 

2.27 

33.1 

3.81 

24.2 

1.516 

1.883 

82.9 

Dec  30 

3.72 

29.9 

5.08 

27.0 

1.547 

1.965 

97.0 

Dec  30 

0.25 

30.7 

5.08 

27.0 

1.547 

1.965 

97.0 

*  Soil  totally  thawed. 
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Table  A8.  c'  and  cj>'  calculated  from  CIV  tests  for  each  soil 
condition. 


Date 

c' 

(kPa) 

I*' 

(degrees) 

Thaw 

depth 

(cm) 

Water 

content 

(%) 

Dry 
density 
( glcm 3) 

Total 

density 

(glcm3) 

Saturation 

(%) 

Nov  9 

16.39 

28.0 

3.81 

20.0 

1.5088 

1.811 

67.9 

Nov  10 

15.36 

20.0 

10.16 

23.4 

1.5456 

1.907 

83.9 

Dec  20 

37.34 

21.6 

60.96 

13.0 

1.6992 

1.920 

58.9 

Dec  27 

46.64 

17.0 

60.96 

9.0 

1.6112 

1.757 

35.6 

Dec  28 

12.88 

23.3 

3.81 

24.2 

1.5152 

1.883 

83.0 

Dec  30 

40.37 

11.0 

8.89 

27.0 

1.5472 

1.965 

97.0 

_ 

Table  A9.  Loading  comparison  for  each  test 
method. 


Normal 

stress 

(kPa) 

Deformation 

rate 

(ends) 

Area  of  applied 
normal  load 
<cn?) 

Triaxial 

34-227 

1.27 

20.6 

Direct  shear 

80-234 

0.08 

31.6 

Shear  annulus 

5-36 

4.2 

126 

Instrumented  vehicle 
wheel 

124-200 

10-40 

387-581 

Figure  A9.  Comparison  of  failure  envelopes  for  each  test 
method.  Solid  lines  mark  the  range  of  normal  stress  for 
each  test. 
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Figure  A1 0.  Calculated  c|>'  with  water  content. 


Figure  All.  Calculated  <|)'  from  the  CIV  test  with 
thaw  depth. 
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Table  A10.  Direct  shear  data  for  Lebanon  Sand  at  constant 
density  of  106  pcf  (cross  section  area  =  4.831  in2;  volume  = 
0.007  ft3),  September  1990. 

Sample  Normal  Peak  shear  Mean  peak 

no.  load  (lb)  W  (%)  load  (lb)  (lb) 

1 

2 

3 

54.7 

Liquid 

44.5 

42.8 

47.0 

4 

81.0 

Liquid 

67.0 

5 

65.0 

6 

69.0 

7 

125.0 

Liquid 

86.2 

8 

80.9 

9 

84.6 

10 

160.0 

Liquid 

113.9 

11 

118.4 

12 

115.8 
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A5.  Chemical  Analysis 


Table  A10.  Chemical 
analysis  for  Lebanon 
Sand,  5  Jan  1997. 


p 

152  Mg/kg 

Na 

0 

K 

10 

Ca 

909 

Mg 

1000+ 

pH 

6.6 

OM 

0.32  % 
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APPENDIX  B.  SOIL  CONDITIONS  FOR  SATURATED  TRIAXIAL 
SHEAR  AND  HYDROSTATIC  TESTS 

The  specimens  for  the  saturated  triaxial  tests  were  prepared  based  on  the  pro¬ 
cedures  outlined  in  Appendix  X,  EM  (Engineer  Manual)  1 110-2-1906  (U.S. 
Army  COE  1970).  These  specimens  were  molded  at  specific  moisture  contents 
and  densities  for  the  drained  and  undrained  tests.  The  specimens  were  then  back¬ 
pressure  saturated,  and  after-saturation  conditions  were  also  determined  for  each 
sample.  Tables  B1  and  B2  are  summaries  of  the  moisture  content  and  dry  densi¬ 
ties  for  each  samples  before  and  after  saturation  conditions  for  both  drained  and 
undrained  triaxial  tests.  The  conditions  for  the  hydrostatic  consolation  test  are 
shown  in  Table  B3.  The  pore  pressure  and  displacement  measurements  were 
recorded  for  the  saturated  undrained  triaxial  test  (Fig.  Bl). 


Table  Bl.  Lebanon  Sand  drained  shear  tests. 


Test 

no. 

Water 
content  (%) 

Dry  density 
(pcf) 

Void  ratio 

2DS1 

Initial  Condition 

9 

108.8 

0.53 

After  Saturation 

18 

108.9 

0.53 

2DS2 

Initial  Condition 

9 

108.7 

0.53 

After  Saturation 

18 

108.7 

0.53 

4DS2 

Initial  Condition 

9 

108.7 

0.53 

After  Saturation 

18 

108.7 

0.53 

6DS1 

Initial  Condition 

9 

108.2 

0.54 

After  Saturation 

18 

110.0 

0.52 

6DS2 

Initial  Condition 

9 

108.3 

0.54 

After  Saturation 

18 

108.5 

0.54 
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Table  B2.  Lebanon  Sand  undrained  shear  tests. 


Test 

no. 

Water 
content  (%) 

Dry  density 
(pcf) 

Void  ratio 

2US1 

Initial  Condition 

13 

109.4 

0.52 

After  Saturation 

17 

109.6 

0.52 

2US3 

Initial  Condition 

13 

109.1 

0.53 

After  Saturation 

20 

110.0 

0.52 

4US1 

Initial  Condition 

13 

114.0 

0.46 

After  Saturation 

18 

115.0 

0.45 

4US2 

Initial  Condition 

9 

108.0 

0.54 

After  Saturation 

20 

110.6 

0.51 

4US3 

Initial  Condition 

9 

107.9 

0.54 

After  Saturation 

20 

109.2 

0.53 

6US2 

Initial  Condition 

9 

107.8 

0.55 

After  Saturation 

20 

109.3 

0.53 

6US3 

Initial  Condition 

9 

108.9 

0.53 

After  Saturation 

20 

111.3 

0.50 

Table  B3.  Lebanon  Sand  hydrostatic  consolation  test. 


Type  of  test 

Hydro 

stress 

(psi) 

Hydro 

stress 

(kPa) 

Vol. 

strain  (%) 

Vol. 

strain 

Coeff.  perm, 
(cm/s) 

Void 
ratio,  e 

Load  (LK3) 

1.2 

8.2 

0 

0.000 

2.7x10"* 

0.52 

Load  (LK6) 

5.6 

38.9 

0.89 

0.009 

2.0x10"* 

0.51 

Load  (LK12) 

11.0 

76.0 

2.20 

0.022 

1.4x10"* 

0.49 

Load  (LK25) 

23.8 

163.9 

3.82 

0.038 

8.9x1 0~5 

0.47 

Load  (LK50) 

53.0 

365.5 

5.39 

0.054 

8.0x1 0^5 

0.45 

Load  (LK10) 

104.4 

720.1 

7.21 

0.072 

7.2x1  O'5 

0.42 

Rebound  (RK50) 

51.8 

356.9 

7.12 

0.071 

8.6x1  O'5 

0.42 

Rebound  (RK12) 

11.9 

82.2 

6.66 

0.067 

8.9x1 0^5 

0.43 

Rebound  (RK3) 

2.6 

17.8 

6.36 

0.064 

7.7x1 0-5 

0.43 
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Figure  B1.  Deviatoric  stress  versus  displacement 
and  pore  water  pressures  for  the  saturated 
undrained  triaxial  tests. 
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Pore  Pressure  (kPa) 

0  10  20  30  40 


Displacement  (mm) 
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Figure  B1  (cont.). 
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Pore  Pressure  (kPa) 
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Pore  Pressure  (kPa) 
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Figure  B1  (cont.) 
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APPENDIX  C.  DATA  FOR  UNSATURATED  TESTS  PERFORMED 
AT  OPTIMUM  MOISTURE  CONTENT  AND  DENSITY 
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Figure  Cl.  Unsaturated  triaxial  tests  (undrained). 
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Table  Cl.  Summary  of  Lebanon  Sand  undrained 
shear  tests  for  samples  at  13%  moisture  content. 


Test 

no. 

Average  cr3 
(confining 
pressure)  (kPa) 

Peak  q 
(CTl  -  o3) 
(kPa) 

Max.  mean 
stress,  p 
(kPa) 

2.1 

14 

57 

33 

2.3 

14 

61 

34 

2.4 

14 

86 

43 

2.5 

14 

101 

47 

4.2 

27 

83 

55 

6.1 

41 

159 

94 

6.2 

42 

121 

82 

Table  C2.  Young’s  Modulus  calculations 


Test  name 

2.4 

2.3 

2.1 

2.5 

4.2 

6.3 

6.2 

6.1 

Confining 
pressure,  psi 

2 

2 

2 

2 

4 

6 

6 

6 

Values  per  load, 
psi 

2766.7 

1940.0 

1660.0 

2200.0 

1780.0 

900.0 

583.3 

1727.3 

1287.5 

1355.5 

1080.0 

1411.1 

1240.0 

541.7 

214.3 

1300 

1250.0 

1111.1 

1150.0 

1125.0 

1277.8 

295.0 

933.3 

750 

1111.1 

1111.1 

1000.0 

1375.0 

Average 

1603.8 

1468.8 

1296.7 

1444.4 

1324.5 

578.9 

576.9 

1259.1 

Average 

1453.4 

1324.5 

805.0 

Total  avg. 

1231.3  psi  (8489.5  kPa) 
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APPENDIX  D.  ABAQUS  CODE  FOR  TRIAXIAL  TEST 
SIMULATION 


**  Simulation  of  triaxial  tests  on  Lebanon  Sand 
Material 

k  k 

**  Sample  size,  2.8"  DIA  x  5.5"  High 

k  ■ k 

*HEADING 

Lebanon  Sand,  Triaxial  test 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

**  Soil  Node  and  Element  Definition 

■ k  - k 


* PARAMETER 


disp=1.375  finches 
press=2  #psi 
cohes=10000  #Pa 
beta=55.8  tdegrees 
ecc=0 . 45 
e_o=0 

youngs=8.5e6  #Pa 
v=0 . 3  6 

press=press* **6895  #Pa 
disp=disp/2* . 0254  #m 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


*NODE ,  NSET=CYL 

1,  0.,  0. 

2,  0.0356,  0. 

3,  0.,  -0.1397 

4,  0.0356,  -0.1397 


*ELEMENT ,  TYPE=CAX4 ,  ELSET=CYLINDER 
1,  1,  3,  4,  2 

■krkrkrkrkrki^^rk-krkrkrkrki^i^rkrkrkrk-krki^i^rkrkrkrkrkrki^^rkrkrkrkrkrk^i^rkrki^rkrkrki^^rkrkrkrk-k-k 


- k  -k 

**  Soil  Material  Definition 

k  k 

*SOLID  SECTION,  ELSET=CYLINDER,  MATERIAL=dpsoil 

* MATERIAL ,  NAME=BAS_EL 

*DENSITY 

1890. 
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*ELASTIC 
<youngs>, <v> 

-k  -k 


*MATERIAL,NAME=dpsoil 

*ELASTIC 


<youngs>, <v> 


*CAP  PLASTICITY 

**  COHESN,  BETA,  BIGR,  EVOLPI,  ALPHA,  BIGK 
<cohes>, <beta>, <ecc>, <e_o>, , 1 . 

*CAP  HARDENING 

**  PB  is  hydrostatic  stress  in  pascals,  EVOLP 
strain  and  unitless** 


8. 

2e 

3,  0 

38 

.  9 

e3 , 

009 

76 

.0 

e3 , 

022 

163 

.  9 

e3 , 

038 

365 

.5 

e3 , 

054 

720 

.  1 

e3 , 

072 

*DENS 

ITY 

1890. 


is  vol. 


'k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k^'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'krk-k-k'k-k-k'k 


**  Boundary  Conditions 

■ k  ■ k 


*BOUNDARY,  OP=NEW 

1,  lrr  0. 

3,  1,,  0. 

3,  2,,  0. 

4,  2,,  0. 


*INITIAL  CONDITIONS,  TYPE=STRESS ,  GEOSTATIC 
CYLINDER, -2  64 . , -0 . 1397,  0 . ,  0 . ,  1 . ,  1  . 

'k'k'k'k-k-k-k'k'k'k'k'k-k-k'k'k'k'k'k'k-k-k-k'k'k'k'k'k-k-k'k'k'k'k'k'k-k-k-k'k'k'k'k'k-k-k'k-k'k'k'k'k-k'k 

**  Step  1,  Geostatic,  LoadCase,  Default 

•k  • k 


*STEP,  INC=500,  UNSYMM=YES 
**,  NLGEOM 

Apply  Hydrostatic  Compression 
* STATIC 

.5,  1.0,  l.E-8, 

*DLOAD,  OP=NEW 
CYLINDER,  P3,  <press> 
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CYLINDER,  P4 ,  <press> 

■k  rk 


*OUTPUT ,  FIELD,  FREQUENCY=1 
*ELEMENT  OUTPUT 
S,E, PEQC 
*NODE  OUTPUT 


U 

*OUTPUT ,  HISTORY,  FREQUENCY=1 
*ELEMENT  OUTPUT,  ELSET=CYLINDER 
S,E, PEQC 

*NODE  OUTPUT,  NSET=CYL 
U,  RF 

*END  STEP 


-k-k'k'k-k-k-k'k-k-k'k'k-k-k'k^-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k-k'k-k-k'k'k-k-k^'k-k-k'k'k-k'k 


**  Load  1.375  inch  displacement 
*  * 


*STEP,  INC=10000 ,  UNSYMM=YES 

Load  to  1.375  in  displacement 

**,  NLGEOM 

** STATIC,  DIRECT 

**0.01 

*  STATIC 

0.01,1.0, .00000001 

*BOUNDARY ,  OP=MOD 

1 ,  2 , ,  -<disp> 

2 ,  2 , ,  -<disp> 

*END  STEP 
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APPENDIX  E.  RESULTS  FOR  TRIAXIAL  TEST  SIMULATIONS 

Initially,  small  cohesions  such  as  3000,  1200,  and  1 .0x10"'  were  used  with 
alternating  beta  values  of  34.4  and  26.6°.  The  graphs  plotted  from  these  data 
showed  that  the  results  grossly  underestimated  the  laboratory  test  data.  The  next 
logical  step  was  to  increase  the  cohesion  while  keeping  the  beta  the  same  as  in 
previous  runs.  Varying  the  cohesion  between  15  and  40  kPa  raised  the  simulated 
deviatoric  stress  to  a  level  that  was  sufficient  to  fit  the  control  data.  However,  the 
slope  of  the  simulated  data  was  very  steep  and  then  very  flat  with  an  abrupt 
maximum,  while  the  slope  of  the  laboratory  data  changed  more  gradually.  To 
correct  this,  the  eccentricity  was  changed  from  a  standard  of  0.45  in  the  previous 
runs  to  0.1  and  0.3  (Fig.  El). 


Axial  Strain 


Figure  El.  Comparison  between  results  from  simulation  runs 
using  higher  cohesions  (15,  20,  and  35  kPa)  with  changed 
eccentricity  (0.1)  and  laboratory  data  for  14  kPa  (2  psi)  confining 
pressure.  The  plots  for  the  28-  and  -40  kPa  (4  and  6psi)  confin¬ 
ing  pressures  exhibited  the  same  behavior  as  the  14-kPa  plot. 

The  data  plots  showed  that  this  was  not  the  answer  to  the  slope  problem;  if 
anything,  this  change  made  the  graphs  steeper  than  before.  In  two  simulation 
tests  the  initial  volumetric  plastic  strain  was  also  changed  to  0.007  from  the  stan¬ 
dard  0.001,  and  again  it  made  no  difference  in  the  plots  (Fig.  E2). 

Finally  the  switch  was  made  from  a  Young’s  modulus  of  1.73e6  to  8.50e6 
(the  final  value  averaged  in  Table  C2).  In  addition,  a  new  beta  value  was  recal¬ 
culated  (55.8°)  based  on  the  ABAQUS  definition  of  deviator  stress,  q.  The 
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eccentricity  was  also  returned  to  its  original  0.45,  and  the  initial  strain  was  set  at 
0.  Six  simulations  were  run  using  these  values  and  a  cohesion  of  1.0x1 0"5,  and 
the  resulting  graphs  fit  almost  perfectly  with  our  control  test  data  (Fig.  E3). 


Axial  Strain 

Figure  E2.  Comparison  of  changed  initial  volumetric  plastic  strain 
to  laboratory  data  and  unchanged  volumetric  plastic  strains  for  28 
kPa  (4  psi)  confining  pressure.  There  were  no  simulations  per¬ 
formed  for  the  change  in  initial  volumetric  plastic  strain  for  14  and 
40  kPa  (2  and  6  psi)  confining  pressures. 


140 


0.00  0.05  0.10  0.15  0.20  0.25  0.30 


Axial  Strain 


Figure  E3.  Comparison  of  new  cohesion  and  beta  graphs  with  labora¬ 
tory  data  for  2,  4,  and  6  psi  confining  pressure  (d  =  I.OxlO-5,  beta  = 
55.8°,  evol  =  0,  eccentricity  =  0.45). 
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However,  when  this  set  of  parameters  was  put  into  the  3-D  simulation,  the 
material  model  faltered  and  would  not  work  correctly.  The  top  layer  of  the  soil 
expanded  erroneously  because  there  was  insufficient  material  cohesion  to  support 
a  tensile  load  (Fig.  E4).  The  final  adjustment  was  to  increase  the  cohesion  to  val¬ 
ues  of  5  and  10  kPa  (Fig.  E5).  While  some  sacrifice  was  made  with  the  fit  of  the 
model  to  the  laboratory  data,  the  3-D  model  ran  smoothly  and  gave  a  satisfactory 
result.  A  final  set  of  parameters  was  obtained  that  gave  a  reasonable  fit  to  both 
lab  and  field  observations. 


The  results  of  the  finite  element  simulations  of  the  triaxial  tests  are  given  in 
Table  El. 


Figure  E4.  3-D  ABAQUS  simulation  of  a  tire  in  soft  soil  with  large 
expansion  caused  by  shearing  of  zero  cohesion  material  (vertical 
exaggeration  for  illustration). 


Figure  E5.  Comparison  of  final  model  parameters  with  laboratory  data  for  2,  4,  and  6  psi 
confining  pressure. 
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Table  El.  Summary  of  triaxial  test  simulations  and  test  results. 


d  (Pa) 

Beta  (°) 

R 

Confining 
press  (psi) 

Young's 

modulus 

Volumetric 
plastic  strain 

Peak  deviatoric 
stress  (kPa) 

3000 

34.4 

0.45 

2 

1.73x106 

0.001 

12 

3000 

26.6 

0.45 

2 

1.73x106 

0.001 

7 

1200 

34.4 

0.45 

2 

1.73x10s 

0.001 

13 

1200 

26.6 

0.45 

2 

1.73x10s 

0.001 

9 

1200 

14.4 

0.45 

2 

1.73x10s 

0.001 

5 

1200 

55.8 

0.45 

2 

1.73x10s 

0.001 

40 

I.OOxlO’5 

14.4 

0.45 

2 

1.73x10s 

0.001 

4 

I.OOxlO-5 

34.4 

0.45 

2 

1.73x10s 

0.001 

12 

I.OOxlO-5 

26.6 

0.45 

2 

1.73x10s 

0.001 

7 

5000 

55.8 

0.45 

2 

1.73x10s 

0.001 

47 

1.00x1  O’5 

14.4 

0.45 

4 

1.73x10s 

0.001 

7 

1 .00x1  o-5 

26.6 

0.45 

4 

1.73x10s 

0.001 

15 

I.OOxlO’5 

34.4 

0.45 

4 

1.73x10s 

0.001 

24 

1200 

34.4 

0.45 

4 

1.73x10s 

0.001 

25 

3000 

26.6 

0.45 

4 

1.73x10s 

0.001 

15 

3000 

34.4 

0.45 

4 

1.73x10s 

0.001 

24 

15000 

26.6 

0.1 

2 

1.73x10s 

0.001 

26 

15000 

26.6 

0.3 

2 

1.73x10s 

0.001 

26 

15000 

34.4 

0.1 

2 

1.73x10s 

0.001 

32 

15000 

34.4 

0.3 

2 

1.73x10s 

0.001 

32 

15000 

26.6 

0.1 

4 

1.73x10s 

0.001 

34 

15000 

26.6 

0.3 

4 

1.73x10s 

0.001 

34 

15000 

34.4 

0.1 

4 

1.73x10s 

0.001 

44 

15000 

34.4 

0.3 

4 

1.73x10s 

0.001 

44 

40000 

15 

0.1 

4 

1.73x10s 

0.001 

51 

15000 

26.6 

0.1 

6 

1.73x10s 

0.001 

42 

15000 

26.6 

0.3 

6 

1.73x10s 

0.001 

42 

15000 

34.4 

0.1 

6 

1.73x10s 

0.001 

55 

15000 

34.4 

0.3 

6 

1.73x10s 

0.001 

55 

7540 

34.4 

0.45 

2 

1.73x10s 

0.001 

22 

7540 

26.6 

0.45 

2 

1.73x10s 

0.001 

17 

3040 

34.4 

0.45 

2 

1.73x10s 

0.001 

16 

3040 

26.6 

0.45 

2 

1.73x10s 

0.001 

12 

15000 

34.4 

0.45 

2 

1.73x10s 

0.001 

32 

35000 

34.4 

0.1 

2 

1.73x10s 

0.001 

57 
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40000 

34.4 

0.1 

2 

1.73x106 

0.001 

64 

20000 

34.4 

0.1 

2 

1.73x10s 

0.001 

38 

15000 

34.4 

0.45 

4 

1.73x10s 

0.001 

44 

7540 

26.6 

0.45 

4 

1.73x10s 

0.001 

25 

7540 

34.4 

0.45 

4 

1.73x10s 

0.001 

34 

30000 

34.4 

0.1 

4 

1.73x10s 

0.007 

44 

45000 

34.4 

0.1 

4 

1.73x10s 

0.001 

82 

40000 

34.4 

0.1 

4 

1.73x10s 

0.001 

75 

35000 

34.4 

0.1 

4 

1.73x10s 

0.001 

69 

35000 

34.4 

0.1 

4 

1.73x10s 

0.007 

69 

15000 

34.4 

0.45 

6 

1.73x10s 

0.001 

55 

7540 

26.6 

0.45 

6 

1.73x10s 

0.001 

33 

7540 

34.4 

0.45 

6 

1.73x10s 

0.001 

45 

35000 

34.4 

0.1 

6 

1.73x10s 

0.001 

81 

40000 

34.4 

0.1 

6 

1.73x10s 

0.001 

90 

I.OOxlO-5 

55.8 

0.45 

2 

8.50x10s 

0 

40 

I.OOxlO-5 

55.8 

0.45 

4 

8.50x10s 

0 

79 

I.OOxlO-5 

55.8 

0.45 

6 

8.50x10s 

0 

120 

I.OOxlO-5 

55.8 

0.45 

2 

8.50x10s 

0.01 

40 

I.OOxlO-5 

55.8 

0.45 

4 

8.50x10s 

0.01 

79 

I.OOxlO-5 

55.8 

0.45 

6 

8.50x10s 

0.01 

120 

1000 

55.8 

0.45 

2 

8.50x10s 

0 

41 

1000 

55.8 

0.45 

4 

8.50x10s 

0 

81 

1000 

55.8 

0.45 

6 

8.50x10s 

0 

121 

5000 

55.8 

0.45 

2 

8.50x10s 

0 

50 

5000 

55.8 

0.45 

4 

8.50x10s 

0 

90 

5000 

55.8 

0.45 

6 

8.50x10s 

0 

130 

10000 

55.8 

0.45 

2 

8.50x10s 

0 

59 

10000 

55.8 

0.45 

4 

8.50x10s 

0 

99 

10000 

55.8 

0.45 

6 

8.50x10s 

0 

138 

Table  E2.  Summary  of  Lebanon  Sand  drained  shear  data. 


Test 

no. 

Average 
ct3  (kPa) 

Peak 

deviatoric 
stress,  q  (kPa) 

Corresponding 

strain 

Max.  mean 
stress,  p 
(kPa) 

Initial 

void 

ratio 

Max. 

volumetric 

strain 

2DS1 

14 

42 

0.18 

28 

0.53 

0.029 

2DS2 

15 

41 

0.24 

29 

0.54 

0.028 

4DS2 

28 

76 

0.16 

52 

0.53 

0.032 

6DS1 

43 

133 

0.19 

88 

0.54 

0.038 

6DS2 

39 

107 

0.20 

76 

0.54 

0.039 
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APPENDIX  F.  MATERIAL  PARAMETERS  FOR  THAWING 
LEBANON  SAND  AND  OTHER  MATERIALS 


The  model  parameters  for  the  thawing  Lebanon  Sand  are  listed  in  Table  FI, 
along  with  model  parameters  from  other  studies  for  comparison. 


Table  FI.  Material  constants  for  Capped  Drucker-Prager  model. 


Material 

parameter 

Lebanon  Sand 

Saturated  Saturated 
drained  undrained 
shear  shear 

1 3%  water 
content 

A-4 

Saturated 

drained 

shear 

Snow* 

(SNOW2) 

McCormick 
Ranch  Sand 

Saturated 
drained 
shear  ** 

Lebanon 

Sand 

Final 

£,  Young’s 
modulus,  kPa 

17250 

17440 

8490 

16000 

1379 

689475 

8500 

v,  Poisson’s 
ratio 

0.32 

0.5 

0.32 

0.33 

0.3 

0.25 

0.32 

(3,  Drucker- 
Prager  angle 
of  friction  (°) 

55.6 

56.4 

54.8 

53 

22.5 

14.6 

55.8 

d,  Drucker- 
Prager 
material 
cohesion,  kPa 

0 

0 

16.7 

16.6 

15 

1190 

10 

R,  cap 
eccentricity 

0.45  to 
0.47 

0.45 

0.45 

1.1*10^ 

0.1 

0.45 

epl 

vol  o 

initial  value  of 
volumetric 
plastic  strain 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0 

K,  flow  stress 
ratio 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

*  Shoop  et  al.  (1 999)  and  Shoop  (2001 ). 

“  HKS  (1996)  and  DiMaggio  and  Sandler  (1971,  1976). 
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The  soil  modeled  is  a  frost-susceptible  fine  sand,  which  was  used  in  full-scale  tests  of  paved  and  unpaved  road  sections  in  CRREL’s 
Frost  Effects  Research  Facility  (FERF).  The  soil  was  subjected  to  a  full  suite  of  saturated  and  unsaturated  triaxial  testing,  using  density, 
moisture,  and  loading  conditions  duplicating  those  experienced  during  the  freeze-thaw  testing  in  the  FERF.  Material  parameters  were 
generated  for  a  capped  Dmcker-Prager  plasticity  model.  These  were  calibrated  in  triaxial  test  simulations  using  the  commercial  finite 
element  code  ABAQUS.  The  material  model  was  then  implemented  in  several  three-dimensional  finite  element  simulations  for 
validation  and  robustness.  The  model  for  Lebanon  Sand  was  compared  to  the  same  model  for  other  granular  materials. 


15.  SUBJECT  TERMS 

Frost  action 

Lebanon  Sand 

Material  properties 
Mobility 

Soil _ 

Spring  thaw 

Thaw  weakening 
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